Electronic spectrum in cuprates within the p-d Hubbard model 
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A microscopic theory for electronic spectrum of the Cu02 plane within an effective p-d Hubbard 
model is proposed. Dyson equation for the one-electron Green function in terms of the Hubbard 
operators is derived which is solved self-consistently for the self-energy evaluated in the noncrossing 
approximation. Electron scattering on spin fluctuations induced by kinematical interaction is de- 
scribed by a model dynamical spin susceptibility. Doping and temperature dependence of electron 
dispersions, spectral functions, the Fermi surface and the coupling constant A are studied. 



Recent ARPES studies revealed a complicated char- 
acter of electronic structure in the copper oxide super- 
conductors which are believed to be caused by strong 
electron correlations in cuprates, as was originally sug- 
gested by Andersoni. Below we report the results of 
electronic spectrum calculations for an effective Hub- 
bard model reduced from the p-d model for the Cu02 
plane in cuprates. In these studies for the first time 
we go beyond the mean-field approximation^ or pertur- 
bation approaclA We have solved the Dyson equation 
self-consistently for the thermodynamic Green functions 
(GFs) and the self-energy derived in the noncrossing ap- 
proximation (NCA), as has been done by us for the t-J 
model^. 



I. MODEL AND DYSON EQUATION 

We consider an effective p-d Hubbard model for one- 
hole states with energy Si = — /i and two-hole p-d 
ginglet states with energy 62 = 2ei + Ueff where ^ is 
the chemical potential and an effective Coulomb energy 
Vef f = Apd = Cp — €d (for details see^) : 
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where X"™ — \in){im\ are the Hubbard operators (HOs) 
for 4 states n, m = |0), |cr), |2) = | ti), cr = ±1/2, a = 
—a. The dispersion of holes is determined by the hopping 
parameters: Uj = tSj,i±a^^^ 4- i' (5j_j±a,±a„ {a-x/y = a 
'- lattice constants). We take Apd ~ 8t ~ 3.2 eV and 
= -0.3t < 0. 

By applying the Mori-type projection technique for the 
GFs Gya(<-t') = {{Xia{t)\xj^{t'))) in terms of the two- 
component HOs {Xj^ = {Xf", Xf°}) an exact Dyson 
equation was derived as described in^ with a self-energy 
(SE) as a many-particle GF. By using NCA for the SE, a 
closed system of equations was obtained for the the GFs 
and the SE: 



where £1(2) (q) are spectra for two bands given by the 
matrix = {{[X,^, H], x]^}) x {{X,,,Xl})-\ The 
SE is the same for the both Hubbard bands: 
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X Imx5c(k-q,i^) Im {Gi(q, z) + G2(q, z)} . 
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Gi{2) (q, tj) = - £1(2) (q) - E(q, u))' 



(2) 



Here the interaction t(q) is determined by the hopping 
parameter iV(z, J/) = (l/2)(tanh(z/2T)+coth(i//2r)) 
and XsM.y) = |((iVq|iV_q))^ + ((Sq|S_q))^ is thc 
charge-spin susceptibility where A'q is the number and 
Sq is the spin operators. 



II. RESULTS AND DISCUSSIONS 

The system of equations (01, ® was solved numer- 
ically for various hole concentrations n = 1 + 5 = 
2{X1" -I- Xf') by using the Matsubara frequency rep- 
resentation at temperature T ~ O.OSi ~ 140 K. Ne- 
glecting charge fluctuations, the spin susceptibility was 
described by the model: Imx5(q, i^) — Xo/[l + ^^(1 + 
7(q))] tanh(i'/2r)/[l 4- (i^/ws)^] where ^ is an anti- 
ferromagnetic (AF) correlation length (in units of a), 
iOs — J = 0.4< is spin-fluctuation energy, and 7(q) = 
(l/2)(cosq2: + cos^y). The constant xo = [3(1 — 
|(5|)/2t^,]{(l/Af)^q[l+C2[l+7(q)]-i}-i is deflnedfrom 

the equation (S,S,} = (3/4)(l - |(5|). 

The dispersion curves given by maxima of spectral 
functions A(k,w) = Bi(k) Ai(k, w) -|- S2(k) A2(k, w), 
where _Bi.2(k) are the weights of the bands and 
Ai(2)(k,Li;) = — (l/7r)ImG'i(2)(k, w), were calculated for 
hole doping 5 = 0.1 — 0.3. The dispersion curves and 
the spectral function for 5 = 0.1 (^ = 2.5) reveal a rather 
flat hole-doped band at the Fermi energy (EE) {lo = 0) 
(Fig.HI Fig. 13 . 

At high temperature T = 0.3t for 6 = 0.1 (^ = 1.0) 
(Fig. 13 Fig. 2J) or in the overdoped region {6 = 0.3) the 
dispersion becomes much larger which proves a strong 
influence of AF spin-fluctuations on the 

With doping, the density of states (DOS) shows a 
weight transfer from the upper to the lower band as 
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FIG. 1: Dispersion curves in units of t along the symmetry ° o.o . 

directions r(0, 0) M{-k, tt) X{it, 0) r(0, 0) for 5 = 0.1 "^'^^ 
and temperature T = O.OSt. 

FIG. 4: The same as Fig.Hbut for T = 0.3t. 




FIG. 2: Spectral function in units of t along the symmetry 
directions r(0, 0) M{tv, tt) X{tv, 0) ~* r(0, 0) for 5 = 0.1 
and temperature T = 0.03t. 



shown in Fig. left panel. The self-energy E(k, re- 
veals an appreciable variation with k and doping close 
to the Fermi level. Figure shows a change of disper- 
sion (kink) in the M X direction at the Fermi level 
crossing for S ~ 0.1. For the coupling constant we get an 
estimation A = wp/wo - 1 ^ 2.4 (A ~ 0.7 for 6 = 0.3). 





FIG. 5: Electronic density of states. 
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FIG. 6: A(k,tj) in the M X direction at the Fermi level 
crossing for 5 = 0.1. 



FIG. 3: The same as Fig.Qbut for T = 0.3t. 



The FS changes from a hole arc- type at 5 = 0.1 to an 
electron-like one at S = 0.3 (Fig. |7|- Fig. EJ). 
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FIG. 7: ACk,Lj = 0) on the FS for 5 = 0.1 



FIG. 9: The same as Fig.[7|but for 5 = 0.3. 
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FIG. 8: The same as Fig. □ but for 6 = 0.2. 



To conclude, the microscopic theory based on HO tech- 
nique for the effective p-d model ^ provides an explana- 
tion for doping and temperature dependence of electronic 
spectrum in cuprates which is controlled by the AF spin 
correlations. Superconducting pairing in the model be- 
yond the weak coupling approximation^ will be consid- 
ered elsewhere. 
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